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Baryon-loops vacuum contribution in renormalized models like the Linear sigma model and 
the Walecka model give rise to large unnatural interaction coefficients, indicating that the 
quantum vacuum is not adequately described by long-range degrees of freedom. We extend 
such models into nonrenormalizable class by introducing an ultraviolet cutoff into the model 
definition and treat the Dirac-sea explicitly. In this way, one can avoid unnaturalness. We 
calculate the equation of state for symmetric nuclear matter at zero temperature in a modified 
a — uj model. We show that the strong attraction originating from the Dirac-sea softens the 
nuclear matter equation of state and generates a vacuum with dynamically broken symmetry. 
In this model the vector-meson is important for the description of normal nuclear matter, 
but it obstructs the chiral phase transition. We investigate the chiral phase transition in this 
model by incorporating deconfinement at high density. A first-order quark deconfinement 
is simulated by changing the active degrees of freedom from nucleons to quarks at high 
density. We show that the chiral phase transition is first-order when quark decouples from 
the vector-meson and coincides with the deconfinement critical density. 
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I. INTRODUCTION 

The equation of state for strongly interacting matter at finite density is needed for understanding 
of neutron stars Q] and heavy ion collision phenomenology P|. One of the most interesting open 
questions for such system is the identification of the appropriate degrees of freedom to describe 
the different phases of matter. One expects that effective degrees of freedom change from baryons 
at low density to quarks and gluons at high density. Therefore, at high density nuclear matter 
undergoes a phase transition to deconfined quarks and gluons and chiral symmetry is restored. Our 
current knowledge about this transition is still very rudimentary since QCD lattice computation 
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cannot yet be performed at large density. In this paper we attempt to describe the first order 
liquid-gas phase transition of nuclear matter, the chiral phase transition, and the transition to 
quark matter in a unified framework. For this purpose, we work within an effective chiral linear 
sigma model coupled to nucleons and quarks. 

The chiral linear sigma model and its extension have been successful in explaining the low 
energy physics of mesons and nucleons in vacuum due to utilizing the concept of spontaneous 
chiral symmetry breaking. However this model has been less successful at finite density. Very 
early studies by Kerman and Miller P| showed that the linear sigma model (LuN) model with 
nucleonic degrees of freedom in fact does not lead to saturating nuclear matter in mean field 
approximation. This failure is due to the fact that in chiral models with Mexican hat potential, 
the medium effect moves the vacuum effective potential towards smaller values. This implies a 
smaller curvature and consequently the mass of the scalar field decreases. The scalar field in the 
linear sigma model plays the role of the chiral partner of the pion and the mediator of medium- 
range nucleon-nucleon attraction. Therefore, as its mass becomes smaller the attraction between 
nucleons becomes stronger. This simple effect may destroy the stability of nuclear matter^. 

Boguta ^ showed, however, that a saturating normal ground state in the LuN model can be 
reproduced by introducing a vector-scalar coupling in a chirally invariant way and generating a 
mass to the vector-meson via spontaneous symmetry breaking. In this approach the saturation 
stems from the cancellation of large repulsive vector mesons and attractive scalar mesons which 
is generated through spontaneous symmetry breaking. However, unacceptably large compression 
modulus of nuclear matter (about K ~ 650 MeV) still remain difficult to overcome within this 
approach. Another short-coming of this approach is the absence of the chiral phase transition, 
since the effective nucleon mass tends to grow at high density P|. There has been some attempt 
to overcome these difficulties by introducing an additional field, the dilaton which is devised to 


effectively incorporate the trace anomaly in QCD j^. However, it was shown that a chirally restored 
phase in the presence of vector-meson coupling to the dilaton is only possible at unphysically high 
values of compressibility {K > 1400 MeV) This paper is an attempt to resolve both problems, 
namely high compressibility and the absence of abnormal phase. We will show that inclusion of the 


^ Note that assignment of a single meson serving in both roles does not nesscesaily lead to a consistent picture. For 
example, in the successful non-chiral Walecka model Hi, the scalar meson only plays the role of the mediator 
of mid-range attraction between nucleons. One should also note that the chiral symmetry realization of the linear 
sigma model is neither unique nor essential. For example, in a non-linear realization |^, the pions field appears 
alone and in vector manifestation scenario of Harada and Yamawaki i, the rho meson is taken as chiral partner 
of the pions. 
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Dirac-sea in a non-renormaliazble fashion reduces the compressibility. The chiral phase transition is 
expected to occur at high density where the baryonic degrees of freedom are no longer appropriate. 
We show that the chiral restoration can be described by incorporating the quark degrees of freedom 
into the LcrN model at high density. We simulate the quark deconfinement at finite density by 
changing the active degrees of freedom from nucleons to quarks at high density. This leads to a 
first order deconfinement phase transition, in accordance with indications from earlier studies [l(] |. 
In our simple model quark deconfinement and chiral restoration are interconnected. A general 
outcome of our analysis concerns the crucial importance of deconfinement for understanding the 
chiral restoration at finite density. 

Initially, the quantum hydrodynamical model was based on renormalizable field theory. There¬ 
fore, a systematic inclusion of vacuum loops was necessary in order to accommodate the response of 
a filled Dirac-sea to the presence of valence nucleons. By renormalizability requirement one could 
then absorb the vacuum loops and redefine the coefficients of the effective potential. However it 
was shown that the resulting effective potential does not provide an acceptable description of the 
properties of finite nuclei, and inclusion of next-to-leading order loop expansion even worsens the 
situation Furnstahl et. al. Q showed also that one-baryon-loop vacuum contribution in 


renormalized models like the LuN model and the Walecka model give rise to large unnatural co¬ 
efficients based on “naive dimensional analysis” proposed by Georgi and Manohar indicating 
that the quantum vacuum is not described adequately by long-range degree of freedoms. 

One possible resolution to the unnaturalness has been non-renormalizable effective held theory 
approach, where one includes all possible terms consistent with the underlying symmetry, hoping 
that short-range physics can be taken into account by non-linear higher order interaction terms. In 
this way, there is no reference to a Dirac sea of nucleons. In this paper, we take a radical approach 
and introduce an ultraviolet cutoff A into the dehnition of the chiral sigma model. In our approach 
the cutoff is taken as new extra parameter which needs to be hxed by some phenomenological input. 
Therefore, our model by default becomes non-renormalizable, similar to the Nambu-Jona-Lasinio 


Q. 


(NJL) model [1^. In this way by construction one does not need to absorb the vacuum loops into 
interaction terms, and therefore the unnaturalness problem can be avoided. With the cutoff we can 
discard the short-range physics which cannot be described by long-range degrees of freedom. In 
other words, the sharp cutoff A separates the active negative energy states from the ones that does 
not contribute. In this way, we explicitly retain the intuitive concept of Dirac sea as a counterpart 
of Fermi sea. 

In the first part of the paper (section II) we study nuclear matter properties within an extended 
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L(tN model in present of explicit Dirac-sea. In the second part of paper (section III), we incorporate 
the quark degrees of freedom above deconfinement density. Having calibrated the parameters of 
the model to the empirical nuclear matter saturation, we investigate the implication of quark 
deconfinement in the equation of state and in the description of the chiral restoration. Some 
concluding remarks are given in section IV. 


II. NUCLEAR MATTER AND NATURALNESS 


We consider the chiral LitN model with exact global SU (2) x SU (2) symmetry, which contains a 
pseudoscalar coupling between pions and nucleons, plus an auxiliary scalar field a. We also include 
a massive isoscalar vector field u)^, 

9sn 4 “ 4 “ C'cj-k 4 “ C,auJ: ( 1 ) 

where Ca-w and are defined as 

+ d^Tid^TT)+ + Ti‘^)+ +Tf‘f, 

C-.. = + vf^) , (2) 

and the field tensor is defined — dyOJfj,. The nucleon mass Mw and vector-meson mass 

rriv at rest are generated through spontaneous symmetry breaking by a-held 

MN{cr) = 9sncr, m^{a) = g^a. (3) 


For uniform nuclear matter, the ground state is obtained by filling energy level with spin-isospin 
degeneracy 7 ^ = 4 up to the Fermi momentum kp, where kp related to the baryon density p by 


P = 


In 


rkp 


(Fk. 


( 4 ) 


(2vr)3 Jo 

We employ the conventional mean-field approximation by expanding the meson fields around their 
expectation values a = {a) and ui^ = neglecting meson fluctuations. The mean values of 
the pion field and space component of the vector field vanish for a baryonic matter at rest due to 
symmetry. Having used the equation of motion for the field, the energy density can be expressed 


as 
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where nucleon and vector masses are defined in Eq. ©• The (To denotes the mean-value of sigma 
field in vacuum. The second line in Eq. © describes the negative-energy nucleon Dirac-sea con¬ 
tribution ^vac subtracted from its corresponding value at zero density. This term is divergent, we 
regularize it by a sharp ultraviolet cutoff A. The linear sigma model is renormalizable, therefore 
any cutoff dependence should be removed by some kind of renormalization scheme in a such way 
that the chiral symmetry of the model is preserved. The vacuum expression can be expanded 
in powers of M^r/A which creates an infinite series in terms of scalar mean-field (notice that the 
only first four orders are divergent). By employing a renormalization scheme, one can then partly 
absorb these terms into the coefficients of the effective potential Q (note also that this procedure 
is not unique), 

4 \ 

( 6 ) 




In aijgsncry ) , 


(dvr) 




i=0 


where /r denotes the renormalization scale (which is normally chosen to be nucleon mass at zero 
density) and Oj denotes the counterterms. The logarithmic term M^ln^^^, however cannot be 
removed in this fashion. The renormalizability requirement imposes very restrictive conditions on 
these highly nonlinear terms. It has been shown that such a renormalized theory is i n ap parent 
conflict with the naturalness property based on Georgi’s naive dimensional analysis [14 1. This 
indicates that the power counting will not be correct^. 

The basic assumption of naturalness is that once the appropriate dimensional scale has been 
extracted using “naive dimensional analysis” proposed by Georgi and Manohar ^ |, the remaining 
dimensionless coefficients should all remain of order unity. If the naturalness assumption is valid, 
then the effective Lagrangian can be truncated within a reasonable error. Notice that there is no 
general proof of the naturalness property, since we do not know how to derive effective hadronic 
Lagrangian from QCD. Nevertheless, phenomenological studies support the validity of naturalness 
and naive power counting rules. 

The first attempt to make this model compatible with the naturalness property has been to relax 
the renormalizability requirement and include all non-linear terms which are allowed according 
to the underlying symmetry 2. Then, one can fix the parameters of the model with some 
phenomenological input. In this fashion, no reference to Dirac-sea is invoked and one retains only 
the valence nucleon explicitly. All vacuum-loop effects will then be hidden in the new interaction 


Moreover, the resulting renormalized effective model fails to reproduce observed properties of finite nuclei, such as 
spin-orbit splitting, shell structures, charge densities, etc m. 
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terms added to the original model. Here, we follow the same line of idea, namely relaxing the 


renormalizability condition by adding the ultraviolet cutoff A to our model definition, but we keep 


Q 


the concept of the Dirac-sea in a similar way to the non-renormalizable NJL model [1^. Here, the 
Dirac-sea is explicitly kept without invoking any new ad hoc interaction terms by expansion. 

The in-medium scalar mean-field a (and its corresponding value in vacuum ctq) or equiva¬ 
lently, the effective nucleon mass ( and M]\f{ao) at zero density) is determined via the 

self-consistency condition by minimizing the energy-density H(cj) with respect to a, 






TT^ 


kFE{kF) — AE{A) + In 


kF + EOvf) 


= 0 , 


(7) 


where E{kF) and E{A) in the above expression are defined by E{k) = + Tbis non¬ 

linear equation should be solved for every point of density associated with Fermi-momentum kF- 
The sigma mass is given by 

d‘^n 


ml = 




( 8 ) 


where the derivatives are evaluated at the scalar mean-field a solution of the gap equation ((T)). In 
the presence of the nucleon Dirac-sea, one does not need to constraint the form of potential from 
outset with < 0 in order to ensure the spontaneous symmetry breaking. The negative energy 
Dirac-sea contribution can indeed generate a dynamical broken vacuum analogous to the NJL 
model. Due to spontaneous symmetry breaking the scalar mean field acquires a nonzero vacuum 
expectation value which is equal to the pion decay constant /jr. Here, the pions are the massless 
Nambu-Goldstone bosons. 

Our model contains six parameters gsn,gv,rnl,^,A and g^n- In contrast to NJL type models, 
here everything can be determined from the gap equation^ ( 0 . In the NJL model the meson masses 
and pion decay should be solved as well since they are described as collective quark-antiquark 
excitations, while in LcrN model they are represented as dynamical fields. We take the empirical 
values for pion decay = 93 MeV, nucleon mass in vacuum Mj-qiao) = 939 MeV and vector-meson 
mass in vacuum m^{aQ) = 783 MeV. These choices fix g^n = Mj-q /= 10.1 and g^ = = 8.42 

via^ Eq. ( 0 . The vector field coupling g^n is treated as a free parameter and will be adjusted in 
medium in such a way that the equation of state reproduces nuclear-matter saturation properties, 
which we define as a binding energy per nucleon E/A = —15.75 MeV at a density corresponding to 


® Notice that in general the gap equation may have many solutions, one has to find out which solution minimizes 
the effective potential. 

The value of gsn differs from the experimental value gsn = 13.5. This is a typical situation in chiral models with 
gA = ^ where nucleons are taken stuctureless. 
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TABLE I: The parameters A, A and gm for sets Al, A2 and B1 — B3 are given (see the text for details). 
All parameter sets reproduce the empirical saturation point {EbIA = —15.75 MeV, po = 0.148 fm“^). The 
corresponding sigma mass rtia in the vacuum and the resulting nucleon mass M]\j(a)\pg = and the 
compressibility K at saturation density are also given. 


Parameter 

set Al 

set A2 

set B1 

set B2 

set B3 

A (MeV) 

256 

300 

324 

444.2 

483.2 

mliGeV^) 

-0.147 

-0.041 

0.026 

0.544 

0.77 

A 

31.5 

27.4 

25 

7.0 

0 

9vn 

6.69 

6.62 

6.54 

6.23 

6.07 

mcr(ao) (MeV) 

806.5 

805 

810.5 

810 

817 

M((f(MeV) 

757.2 

761.2 

766.3 

781.3 

784.5 

K (MeV) 

490 

478 

455 

396 

370 


a Fermi momentum of kp = 1-3 fm~^. The coupling A and A can be determined by requiring 
the value of the pion decay /tt = 93 MeV and value of sigma meson mass dehned in Eq. 0. In 
this way, one still cannot uniquely fix the value of the cutoff. Therefore, for sake of generality we 
consider various parameter sets Al,A2 and B1 — B3 given in table I which are determined with 
the above-mentioned procedure and all pass through the empirical nuclear matter saturation point 
{EbIA = —15.75 MeV, po = 0.148 fm“^). The empirical compression modulus (or pressure) is not 
used in the fitting procedure. Generally, increasing the cutoff shift the saturation point to higher 
density, however for not very high cutoff (given in table I) one can adjust other parameters in order 
to reproduce the empirical saturation point. Parameter sets Al and A2 are chosen with condition 
that < 0, (i.e. with initial Mexican hat potential even without inclusion of the Dirac-sea) while 
for parameter sets B1-B3 we have > 0 (in this case, the Dirac-sea generates a Mexican hat 
potential). For various parameter sets a sigma-meson of mass about 810 MeV is needed in order 
to reproduce the empirical saturation point (see table 1). The value of nucleon masses at the 
saturation point for various parameter sets are also given in table L 

On the right panel of Fig. 1, we show the mean-value of scalar field a with respect to density 
p/po- As it is observed as the density is increased the mean-value of scalar field decreases. However, 
at high density above 2 po the curves (for different parameter sets) bend upward, indicating that the 
chiral restoration can not be described in this model 0,0]. This is due to the fact that the model is 
inadequate in short-distance physics while the meson degrees of freedom a and co account for nuclear 
interaction at short distance. Moreover, the vector-meson mass appears in the denominator of the 
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FIG. 1: On the right panel, we show in-medium pion decay as function of density p/po (with nuclear matter 
density po = 0.15/m“^) for various parameter sets given in table I. On the left panel, we show the binding 
energy per baryon Eb/A as a function of density p/po- 

energy density Eq. © and cannot vanish. On the other hand, the vector-meson mass is locked 
to the scalar field by chiral symmetry. Therefore, a as well cannot vanish, our numerical results 
reflect this fact. In the next section, we will introduce a scenario in which the chiral restoration 
problem can be resolved. 

In Fig. 1 left panel, we show the density dependence of the energy per baryon EsjA = p — 
Misf{d'Q) calculated for different parameter sets given in table I. The value of O is calculated from 
Eq. 0 by making use of the self-consistency equation 0. It is observed from Fig. 1 that as we 
increase the cutoff (increasing the contribution of the Dirac-sea) the equation of state becomes 
softer. In order to find out quantitatively the stiffness of the equation of state, we compute the 
compression modulus. It is defined as 


K = 9 


j' 


dp 


(9) 


The compression modulus at the saturation density for parameter set B3 and Al are 370 MeV 
and 490 MeV, respectively. The compression modulus of other parameter sets lie between these 
two values (see table I). These values are still bigger than the empirical one K = 200 — 300 MeV 
17l| . Nevertheless, it improves compared to no-sea approximation where we have K ~ 650 MeV 
J Note that it has been shown by Prakash and Ainsworth that the empirical compressibility 
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FIG. 2: On the left, the pressure with respect to the energy is shown for parameter sets Al and B3 given in 
table I. On the right side, we show the pressure as a function of chemical potential /r for the same parameter 
sets. 


can be obtained when the vacuum-loops in a renormalizible fashion are included. However, this 
results to unacceptably large scalar sigma mass of about 1 GeV and unnaturalness of the effective 


Q 


Lagrangian as we already discussed. It has been pointed out by Furnstahl and Serot IlJ that the 
scalar meson mass must be light enough in order to avoid strong fluctuations in the charge density, 
which signal impending instabilities. 

By increasing the cutoff, the value of the nucleon mass at saturation density slightly increases. 


For the various parameter sets given in table I, at saturation we 


lave M'^/Misf{cro) = 0.81 — 0.84 


3 


which is in good agreement with analysis by Mahaux and Sartor 

Next we compute the pressure P at zero temperature as function of density p, the baryonic 
chemical potential p and energy. We use thermodynamical identities 



P = pp — VL = 


n d 

PTp 



( 10 ) 


In Fig. 2 left panel, we show the pressure P as function of energy E = Vl for set Al and H3 (other 
parameter sets behave similarly). The causal limit E = P is also shown in the figure. It is obvious 
that all the parameter sets studied here respect the causal condition dP/dE < 0. Therefore, non- 
renormalizibility in our model does not lead to conflicts with causality. The pressure at vacuum and 
saturation density is zero P{po) = 0 which gives a nontrivial constraint on the model parameters. 
The phase transition becomes more apparent by computing pressure as a function of chemical 
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potential fi. A first order phase transition is manifested by appearance of several branches of Pip)- 
In Fig. 2 right panel, the pressure with respect to chemical potential for parameter sets A1 and 
S3 are shown. These curves show the liquid-gas phase transition in nuclear matter. The first 
order phase transition occurs at Hc — Mn{^o) — 15.75 ~ 923.25 MeV for all various parameter sets 
given in table I. At this point, the slope of curves P{p) exhibits a jump corresponding to difference 
of densities between coexisting phases. At a hxed chemical potential only the highest pressure 
corresponds to a stable phase. 

Some remarks here are in order. One can observe that as we increase the cutoff, i.e increasing 
the negative energy Dirac-sea contribution, a stronger coupling ml (with positive sign) is needed, 
while the coupling A subtantialy decreases (see table I). This can be understood since the Dirac- 
sea contains many-body forces and in principle there is no need to add arbitrary extra many-body 
forces in terms of potential into system when the Dirac-sea is enhanced. On the other hand, 
the nuclear matter saturation properties can be better described (compressibility is reduced) with 
increasing cutoff in a acceptable range. The Dirac-sea also generates a vacuum with dynamically 
broken symmetry. Note that the exact connection between the chiral dynamics and nuclear matter 
saturation mechanism is not yet completely understood. Recently we have shown in a relativistic 
i^addeev approach Q that nuclear matter saturation might be related to partial chiral restoration 

3 - 

In the bosonized version of the NJL model, self-interaction term is absent A = 0 (even 
in vacuum), therefore our cutoff LitN model calibrated to nuclear matter saturation properties 
becomes very similar to the NJL model since we have A —> 0 as we increase the cutoff (see table I). 
It is interesting to note as well that in mean field approximation, the nuclear matter in the NJL 


model does not saturate without inclusion of a vector-scalar fermion self-interaction 


221. Therefore 


the same saturation mechanism is behind both models. Finally, as can be noticed from table I, 
the acceptable range of cutoff is too low A ~ 250 — 485 MeV. A similar r ang e for the cutoff has 
been also reported in the NJL model with nucleonic degrees of freedom 22, [2^ • This is in contrast 
with the case that the model is dehned in terms of quarks where the allowed range of the cutoff is 
typically larger. This may indicate that by changing the relevant degrees of freedom from nucleons 
to quarks, the model can be extended to higher energy. 
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III. CHIRAL RESTORATION AND DECONFINEMENT 

In this section we consider the implication of quark deconfinement in the nuclear matter equation 
of state. We incorporate the deconfinement effect by adding the quark dynamics to the nucleonic 
Lagrangian Cn defined in Eq. o and assume that at a critical Fermi momentum kc (corresponding 
to a critical density pc), the baryonic degrees of freedom are replaced by quark ones, while, keeping 
the description of scalar and pion fields intact at higher density. We cannot derive the exact value 
of deconfinement density pc in the context of mean field theory where the internal structure of 
hadron is neglected, therefore we do not constrain ourself to only one value for pc and consider 
the implication of different values for pc- The energy density can be written 

Q. = VLn + ^qO{p - pc), (11) 

where Q-n is obtained from defined in Eq. o and Hg is obtained from the following quark 
Lagrangian 

C-q = - Qsq (o' + *75^7?)]V'g + C (TTT “1“ (12) 

where and are defined in Eq. ©. We do not couple the vector-meson field directly to the 
quark field since the vector-meson in this model works against chiral restoration phase transition 
at high density. Therefore, we assume that the quark field is only coupled to the scalar cr and the 
pion vf fields in chirally invariant way. But we allow the vector-meson field to couple to the scalar 
field through interaction term Since quarks are coupled to scalar field, the vector field is 

not completely decoupled from quarks. The quark mass like the nucleon and vector meson masses 
Eq. (| 31 ) are generated through spontaneous symmetry breaking via scalar field, 

mq{a) = Qsqcr. (13) 

The energy density at finite density in mean field approximation as a function of Fermi momen¬ 
tum can be obtained as 

^ ~ 9n‘^m^{a) ~ — kc)) + (^^ — d-g) + ^ “ ^o) 
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where the quark mass the nucleon mass and the vector-meson mass rriy are defined in 
Eqs. (Idllhf) . 7q = 12 and 7^ = 4 are the degeneracy factor for quarks and nucleon, respectively. We 
have also included the corresponding Dirac-sea of the nucleons and the quarks as in the previous sec¬ 
tion. A priori we do not know the importance of the vector-scalar interaction above deconfinement 
density. We consider two cases above the critical deconfinement density pc'- Qv = ^ corresponding 
to totally decoupling of the vector field from matter and Qv ^ ^ when the vector-scalar interaction 
is present®. In the energy density equation (d, the parameter TZ can have two values 7^ = 0,1 
corresponding to the solutions when = 0 and Pv ^ 0 above the critical deconfinement density pc, 
respectively. Note that the parameter TZ was merely introduced in order to write both solutions in 
one equation. 

The Fermi momentum kp can be related to baryonic density by 

2 k^ 

P = ^- 

The mean-value of scalar field a is determined by in-medium gap equation: 

~ 97]-4^^^3 — kp) + TZk^ 6 (kp — kc)) + Act® -I- fq{kc) — 

- fn{kc) + fn{kp) 0 {kc - kp) + {fn{kc) + fq{kp) “ fq{kc)) (^{kp - kc) = 0 , (16) 



with 


Uk) = 


ixglx<^ 


ky/k"^ + {gsx(y)‘^ - iPaxCr)"^ In 


k + \/k^TJg^^^ 


(17) 


{gsxcr)"^ J ’ 

where we used the notation that x = q,n stands for the quark and the nucleon terms, respectively. 

The quark interaction term Eq. m introduces a new extra parameter psq which is fixed by 
choosing the quark mass at vacuum mq(cJo) = 300 MeV, therefore we have psq = 3.23. The other 
parameters are determined for a given value of deconfinement density pc (or associated Fermi 
momentum kc) in the same procedure described in section II, so as to reproduce the empirical 
nuclear matter saturation point {EpIA = —15.75 MeV, po = 0.148 fm“®). Therefore, we calibrate 
all parameters to nuclear matter properties. We choose again, Mj\f{ao) = 939 MeV, my = 783 MeV 
and = 93 MeV in vacuum. Notice that, here the cutoff A depends on the deconfinement density 
or equivalently Fermi momentum kc, see Eq. (d. This is in fact a cutoff available for a quark in 
the phase space measured relative to the critical deconfinement Fermi momentum kc- A density 


dependent cutoff was also invoked by Harada, Kim and Rho m to stu^ the chiral restoration in 


the context of vector manifestation scenario of Harada and Yamawaki 


4 


® In this case, the value of the coupling is fixed by the empirical vector-meson mass in vacuum. 
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TABLE II: The parameters A, A and g^n are shown for sets Cl and C2 for two values of the critical 
deconfinement density pc- We assume a cutoff A = 483.2 MeV and the coupling psq = 3.226 (which 
yields constituent quark mass 300 MeV) for all parameter sets. All parameter sets reproduce the empirical 
saturation point (Eb jA = —15.75 MeV, po = 0.148 fm“^). The corresponding sigma mass Too- in the vacuum 
and the resulting in-medium nucleon mass at saturation density are also given. 


Parameter 

set Cl 

set (72 

kc (MeV) 

368.8 {pc = 3po) 

292.7 {pc = 1.5po) 

m2(GeV^) 

0.48 

0.376 

A 

15 

20 

Qvn 

6.35 

6.54 

mcr(CTo) (MeV) 

815 

807 

M*jq{yieY) 

772.2 

762.3 


In contrast to the case that transition to quark degrees of freedom is not incorporated, here the 
cutoff can be taken larger and cutoff-dependence is less pronounced. In order to compare our result 
with the previous section we take the cutoff given in table I for set i?3, A = 483.2 MeV. In table 
n, we show two sets of parameters determined for the fixed cutoff A = 483.2 MeV assuming that 
deconfinement take place at density pc = 3/9o and 1.5po for parameter sets Cl and (72, respectively. 
As we already pointed out the cutoff depends on the deconfinement density, this means for a fixed 
cutoff one needs to adjust other parameters in order to reproduce again given phenomenological 
input, see table II. 

In Fig. 3 left panel, we show the mean-value of scalar field as a function of density for various 
parameter sets S3 and (71, (72 in the presence of the vector-scalar interaction terms 7 ^ 0 given in 
table II. For comparison, we also plotted again the corresponding result for parameter set S3 given 
in table I (without the deconfinement effect). It is seen that exactly at the critical deconfinement 
density pc, the slope of in-medium scalar-mean field changes toward chiral restoration (for set 
(72 this is more obvious). This is in contrast to pure nuclear matter case set S3 (see also Fig. 1 
left panel) where chiral restoration seems to be in conflict with nuclear matter properties. On 
the right panel of Fig. 3, we show the binding energy per baryon for various sets. It is seen that 
the equation of state at very high density p ^ Qpo become softer than pure nuclear matter. At 
the deconfinement density there is a pronounced kink which is more obvious if the deconfinement 
occurs at lower density (see Fig. 3, set (72). This is due to the fact that a sharp boundary is 
assumed between the baryonic and the quark phases. In principle, instead of the theta function in 
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FIG. 3: On the left panel: the scalar mean-field values with respect to density p/po (with nuclear matter 
density po = 0.15 /to“^) is shown when 7 ^ 0 for various parameter sets B?, and (71,(72 given in table I 
and II, respectively. On the right panel: we show the binding energy per baryon Eg jA as a function of 
density p/po for the same sets of parameters. The vertical dotted lines show the positions of the critical 
deconfinement density pc- 


Eq. (11411 one may employ a smooth function. Note also that here the formation of Cooper pairs 


and color-superconductivity is not taken into account 




It has been shown that such effects 


might have important consequences for the transition to quark matter within compact stars 


25|. 


In particular, diquark degrees of freedom seem to be essential for the emergence of coexistence 


Q- 


region between the hadronic and deconfined phases 

Next, we switch off completely the vector-scalar interaction term = 0, (i.e. 7^ = 0 in 

Eq. (I14l) i above the deconfinement density pc- In Fig. 5, we show the mean-value of scalar field as 
a function of density for parameter sets Cl and C 2 when Qv = ^ for p > pc- If the deconfinement 
take place at higher density the chiral phase transition will be stronger and coincide with the 
deconfinement critical density. In this case the chiral phase transition is first order. While for a 
soon deconfinement (for example, parameter set C2) the chiral symmetry is partially restored at 
the deconfinement density and full chiral restoration postpones at relatively higher density than 
the deconfinement density pc- It is interesting to note that the soon deconfinement (parameter set 
Cl) is not compatible with Brown-Rho scaling 27l |. 


f-K{cr) _ _ rny{a) 

U{ao) Mn{(Jo) my{ao)' 


( 18 ) 
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FIG. 4: We show the scalar mean-field values as a function of density p/po for parameter sets Cl and 
C2 given in table II when the vector-coupling is zero above the critical deconfinement density pc- 


since we have my{a) = 0 above the deconfinement when = 0, but as it can be seen from Fig. 5 for 
parameter set Cl, the mean-value of scalar field 7,^ above the deconfinement is partially restored 
and is not zero. If we assume that Brown-Rho scaling is valid at various density, we can put a 
constraint from below on the onset of critical deconfinement density pc in our minimal Lagrangian 
model. We found that for the deconfinement density pc > 2pQ, we have fn{pc) = 0 and my{pc) = 0 
at the same time consistent with Brown-Rho scaling and vector manifestation scenario which is 
considered to be a generic feature of effective field theory matched to QCD y. Therefore, for 
Pc > 2p0) the chiral restoration is first order and coincide with the critical deconfinement density 
PC- 


IV. CONCLUSION 


We investigated the implication of the baryonic Dirac-sea in the presence of an explicit ultravi¬ 
olet cutoff in the nuclear matter equation of state within the LcrN model. In this way, the model 
becomes non-renormalizable and one can avoid the unnaturalness problem [l^ . We investigated the 
parameter space of the model as a function of cutoff which reproduces the empirical nuclear matter 
properties. Note that the empirical compression modulus was not used in the fitting procedure. 
We showed that a filled Dirac sea of baryons can produce a strong attraction in nuclear matter. It 
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softens the equation of state and it generates a vacuum with dynamically broken symmetry even 
without a presence of a negative mass term in the linear sigma model. This also clarified a possible 
connection between the chiral LcrN model and bosonized version of the NJL model in a nuclear 
matter medium. 

Although the LciN model in the presence of the vector-scalar interaction can describe the nuclear 
matter saturation, it is unable to accommodate the chiral restoration phase transition with /^r = 0. 
The role of vector-meson is very subtle in this model. On one hand it is needed in order to reproduce 
the empirical saturation properties, but it works against the chiral restoration at high density. In 
order to study the implication of the first-order deconfinement effect in the equation of state and 
chiral phase transition, we simulated the deconfinement effect by changing the active degrees of 
freedom from nucleons to quarks above the critical deconfinement density. At the moment, we do 
not know yet the order of a possible deconfinement transition and its critical density. However, 
there are some indications that it might be first order and takes place at about several times 


the nuclear matter density [lli| . Here, we considered the implication of different choices of the 
critical deconfinement density. We investigated the role of the vector-meson and the deconfinement 
effect in describing the chiral phase transition when our model parameters are calibrated to the 
nuclear matter properties. The nature of the chiral phase transition is sensitively dependent on the 
presence of the vector-scalar interaction. The vector-meson has been reported by many authors to 
be important for realization of the chiral phase transition Q . We showed that if the vector-meson 
is completely decoupled from the quark matter Qv = ^ above the critical deconfinement density, 
chiral phase transition can take place and its position depends on the given critical deconfinement 


density. We showed in order to be consistent with Brown-Rho scaling m and vector-realization 


3 


based on hidden-local symmetry, the critical deconfinement density must be bigger than 2po 
where pQ is the nuclear matter density. Then the chiral restoration is first order and is coincident 
with the deconfinement density®. 

It is of interest to extend scheme presented in this paper to finite temperature and also study 
the implication of this nuclear-quark equation of state on neutron star. 


® Notice that in contrast with the finite density where Lattice simulation is plagued with sign problem, lattice 
calculation has been very fruitful for finite-temperature (at zero baryon density) studies An analysis of the 
lattice simulation shows that chiral symmetry restoration and deconfinement phase transition occurs at the same 
temperature |23. |3(jl. It has been also shown that the vector coupling is small compared to the scalar coupling at 
high temperature 1,3IT . 
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